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G. Alvin Pierce and Anand R. Vaidyanathan 
Georgia Institute of Technology 


SUMMARY 

A variety of approximate methods are available for the estimation of airloads on 
helicopter rotor blades. These methods vary widely in their degree of approximation, 
accuracy and detail of prediction. One such method has been suggested by Van Holten 
which uses an acceleration potential description of the flow field and a matched asymptotic 
expansion technique to calculate unsteady, three-dimensional airloads on a rotor blade in 
forward flight. 

The study presented here was undertaken to examine the theoretical basis and 
computational feasibility of the Van Holten method, and to evaluate its performance and 
range of validity by comparison with experiment and other approximate methods. It is 
found that, within the restrictions of incompressible, potential flow and the assumption of 
small disturbances, the method does lead to a valid description of the flow. However, due 
to the nature of the basic assumptions, the method begins to break down under conditions 
favoring non-linear effects such as wake distortion and blade/rotor interaction 

INTRODUCTION 

This study is concerned with the calculation of three-dimensional, unsteady airloads 
on a helicopter rotor blade in forward flight. A detailed knowledge of these airloads is of 
importance in problems of noise and vibration reduction and in dynamic and aeroelastic 
analyses of the rotor blade. Due to the complexity of the flow field a completely general 
solution to the problem has not been obtained. There are a large number of aerodynamic 
analyses available in which the actual flow field is simplified considerably to make the 
problem solution feasible. Most of these methods are described in reference 1 and, with 
particular reference to the vortex representation in reference 2. The methods range in 
complexity from simple blade element representations to lifting surface models with freely 
distorted vortex wakes with associated ranges in computational expense, accuracy and 
detail of solution. 

A commonly used method is based on a combination of simple momentum theory and 
blade element description, in which a momentum balance is made at each blade section in 
order to determine an effective incidence, from which force and moment coefficients are 
determined. Induced velocity over the disk is usually assumed constant, equal to the value 
obtained from simple momentum theory. Sometimes linear variations of induced velocity 
across the disk are also assumed^ The method does not incorporate three-dimensional and 
wake effects on the airloads other than by the use of empirical factors. The blade element 
description is often quasi-steady. Unsteady aerodynamics can be included by using 
Theodorsen's results (ref. 3) or, with the effect of forward speed variation, Greenberg's 
results (ref. 4). Real fluid effects are accounted for using wind tunnel test results 
appropriate to the Mach number and Reynolds number of the blade section. In spite of these 
shortcomings, this model is perhaps the most widely used in practice, for its main 
advantages are twofold. 



(1) The blade element representation leaves the model open for incorporation of 
empirical factors associated with various effects such as reversed flow, dynamic stall, 
compressibility, etc. (usualiy determined from two-dimensional wind tunnel tests). 

(2) The relative simplicity of the model makes it ideal for inclusion in a more 
extensive helicopter anaiysis. An example of such a model is that of Gormont (ref. 5). 

In order to account for the rotor wake, many methods modei the blade by a simple 
lifting line of bound circulation and then calculate the non-uniform induced velocity field 
generated by the wake of the lifting line, to various degrees of approximation. Finite-chord 
effects are then represented by a blade element model. In the method due to 
Willmer (ref. 6), the wake immediately behind the blade is replaced by a plane, semi- 
infinite wake similar to that of a fixed wing. The remainder of the wake is represented by a 
series of plane, infinite layers underneath the blade. For simplicity, quasi-steady 
assumptions are made so that the wake consists only of trailing vorticity. The wake can 
also be represented by a rigid, helical surface of shed and trailed vorticity from the 
rotating lifting line. Calcuiation of the velocity induced by such a wake at a point on the 
lifting line invoives numerical integration with a singularity in the integrai, which makes it 
necessary to extract the finite part of the integral by some means. This difficulty arises 
because replacing the biade by a lifting line is not a physicaliy acceptabie assumption in the 
vicinity of the blade itself. One way of overcoming this difficulty is to calculate the 
induced velocity along the three quarter-chord line. The only justification for this 
procedure is the observation from steady, two-dimensional flow that, with the distributed 
vorticity lumped into a bound vortex at the quarter-chord point, the boundary condition is 
satisfied exactly only at the three quarter-chord point. Another method of overcoming the 
singularity is to take the time (azimuth) average of the effect of the instantaneous wake 
configuration (ref. 2). This leads to a simplier model which could also have been arrived at 
by replacing the finite-bladed rotor with an actuator disk. This model consists of an infinite 
number of blades each carrying an infinitesimal load, so that the wake of a finite number of 
helical sheets is replaced by a skewed, semi-infinite cylinder fiiled with vorticity. Another 
method consists of considering the continuous distribution of vorticity in the helical wake 
to be lumped into elements of finite circulation at convenient radial and azimuthal 
intervals (ref. 7). If the bound circulation on the blade is approximated by a stepped 
distribution and if the blade is moved impulsively from one azimuth station to the next, 
then the resulting wake is just a network of straight-line segments, a discretized version of 
the continuous wake. The contribution of each segment is known directly as a function of 
its position; hence the wake integration becomes greatly simplified, since it consists only of 
a finite summation. The model can be improved to better account for unsteady 
aerodynamic effects by retaining a continuous sheet for the "near" shed wake (ref. 8). 

Under conditions of low inflow through the rotor and low forward speed, the concept 
of a rigid wake, (i.e., one in which all elements of the wake are convected downward at a 
constant speed to give the wake a rigid helical shape, is not an acceptable one. Interaction 
between the blade and the wake, self-induced distortions in the wake and other nonlinear 
effects become significant. In addition, there is always the process of rapid roll-up of the 
vortex sheet near the blade tip and the consequent effect of a strong tip vortex on the 
same blade and the following one. Miller (ref. 9) introduces the concept of a "semi-rigid" 
wake, in which the downward velocity with which the wake elements are shed varies with 
azimuth but not with radial position. He divides the wake into a near wake, for which the 
blade has a surface representation, and a far wake, for which the blade is a lifting line. The 
bound circulation distribution is also simplified to a constant so that the trailing vorticity 
consists only of a root and tip vortex. It is also possible to carry out a distorted wake 
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analysis by starting either with a rigid wake or some given initial wake configuration, using 
the non-uniform induced velocity field to define a distorted wake for the next iteration and 
repeating the process until convergence is achieved. An example of such a "free" wake 
analysis is that of reference 10. For such an iterative process, the computation time and 
expense for a single wake calculation must be small and the segmented wake model is best 
suited to this need. 

Another approximate method is the "local momentum" approach of reference 11. The 
blade is treated as a series of elliptical wings, each of which contributes an induced 
velocity that is constant along its span. This is combined with a momentum balance at each 
blade element and an "attenuation coefficient" to account for the timewise change of 
induced velocity at a point after blade passage through it. 

In addition to the vortex methods, there are also techniques based on the use of an 
acceleration potential. Dat (ref. 12) models the blade by a lifting line of acceleration 
potential doublets along the quarter-chord line of the blade and satisfies the normal 
velocity boundary condition along the three quarter-chord line. The details of the 
computational procedure are presented by Costes (ref. 13). Although the acceleration 
potential formulation leads to potential discontinuities being confined to the blade surface, 
the computational expense involved is no less than in the vortex (velocity potential) 
formulation. This approach does have the merit that it can be rigorously extended to 
compressible flow; however, computations (ref. 13) seem to show that results for 
compressible flow are well approximated by scaling incompressible flow results. 

In considering the full range of available methods, it is seen that even the distorted 
wake representations are not acceptable under all circumstances. It is possible to account 
for wake distortion and have good correlation with experiment by making use of 
experimental information to define a "prescribed", distorted wake. However, to precisely 
define the prescribed wake for various flight conditions requires an extensive data base of 
experimental results. This has been achieved so far only for the case of hovering flight 
(ref. 14). To extend the data base to forward flight would require entensive experiments to 
cover possible variations of flight conditions. It has been shown (ref. 15) that the use of a 
prescribed wake does produce acceptable results but the results are highly sensitive to the 
prescribed wake geometry. 

The use of potential vortex filaments to model the rotor wake sometimes leads to 
excessive, unrealistic wake effects in the computation, which can be avoided only by 
accounting for factors such as a finite core radius for the actual vortex tube. In view of the 
number of factors that cannot be rigorously accounted for in a rotor aerodynamic 
representation, most methods in practice rely on empirical input in one form or another 
(ref. 16). Any new method proposed for calculating the airloads on a rotor blade in forward 
flight must therefore have its performance assessed against this background. 

A method has been proposed by Van Holten (ref. 17) for calculating the unsteady, 
three-dimensional airloads on a helicopter rotor blade in forward flight. The method uses 
the acceleration potential formulation, together vdth a matched asymptotic expansion 
technique, to generate a solution accurate to 0(A" ), where A is the aspect ratio of the 
blade. This method has not been studied in sufficient detail to determine its value relative 
to other available techniques. This program has therefore been undertaken to carry out the 
following with respect to the asymptotic method. 

(1) To study the theoretical basis and limits of validity 

(2) To verify the expressions derived for computation 

(3) To write a computer program for the computational scheme which calculates 
the pressure distribution on a rotor blade in forward flight 
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(4) To apply the program to calculate results for some cases for which 
experimental results are available for comparison 

(5) To calculate results for the same cases using other approximate methods so that 
results may be compared relative to computational expense 

(6) To consider the possibility of using measured section properties 

(7) To study the efficiency of the computational scheme 
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Notations: 
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( )* 


cylindrical coordinate 

spanwise coordinate along blade with hub as origin 

semi-span of rotor or fixed wing 

time 

freestream velocity 
flow velocity components 
flow velocity 

flow-axes coordinates (fig. 1) 

fixed-wing coordinates (figs. 2, 3 and 4) 

blade-axes coordinates (figs. 1 and 2) 

rotor-axes coordinates (fig. 1) 

tip path plane incidence, forward tilt positive 

blade flapping angle 

blade inertia coefficient (Lock number) 

linear twist, root pitch angle - tip pitch angle 

elliptic coordinate 

prolate spheroidal coordinate 

blade root pitch angle 

variables of integral (eq. (G7)) 

advance ratio, U/fiR^ 

air density 

elliptic coordinate 

cylindrical coordinate 

prolate spheroidal coordinate 

azimuth angle 

rotor angular velocity 

oscillation frequency 

perturbation quantity 
complex amplitude 
non-dimensional 



THE ASYMPTOTIC METHOD 
Analytical Description 


The asymptotic method as proposed by Van Holten (ref. 17) for the determination of 
unsteady, three-dimensional airloads on a helicopter rotor blade is based on the following 
assumptions: 

(a) Incompressible, potential flow 

(b) Disturbances induced by the rotor are small compared to the forward speed 
Each blade of the rotor is presumed to be rigid and rectangular in planform. 

Coordinate systems which are used to describe the rotor configuration are illustrated in 
figure 1. The "rotor-axes" coordinates, x , y , z , are oriented such that the x -y plane is 
the tip-path plane and the freestream direction is parallel to the x -2 plane."^ Tfie "flow- 
axes", X, y , 2, are inclined to the rotor-axes by the angle a > which is £ rotation about y , 
such that the freestream velocity is in the negative x direction. A third system which fs 
fixed to the blade has its origin at the quarter-chord of the mid-span. These "blade-axes", 

^Jb’ ^b’ ^b’ inclined to the tip-path plane by the coning angle such that z. coincides with 
the quarter-chord line. ° 

The momentum equation can be written as 

^ = grad(-p/p) 

Introducting perturbation quantities V' and p' so that V = U + V' and p = Pj„ + p' and using the 
assumption of small perturbations, it can be shown that the p>erturbation° pressure satisfies 
the Laplace equation. 
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The complete linearized boundary value problem for the pressure field around the rotor can 
now be stated: 


2 2 
2 2 

3x‘‘ 3y^ 



( 1 ) 


On the blade surface, the pressure gradient in the 2 direction must attain the value 
(Appendix A) 


such that 
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2 2 2 
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and 

p'-»- - 00 along the leading edge such that the velocity component in .the z direction 
along the quarter-chord line is 
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( — - ) 

i2Rj \ = 0 


(3) 
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The solution to this problem will be divided into a singular part which satisfies = 0 
(singular at the leading edge) and a regular part which is non-singular and satisues the 
given pressure gradient condition of equation (2). 

Singular solution. - In the near field (vicinity of the blade) it is assumed that spanwise 
variations are characterized by a length of the order of the span while the characteristic 
length for variations in the other two directions is of the order of the chord. The governing 
equation for the near field can be written in terms of the nondimensional coordinates 

==5 = ==b = V® 


as 


ajp!. 3 V 1 3 V 

*2 *2 2 *2 ” ^ 


Assuming the near field solution to be described by an asymptotic expansion in inverse 
powers of the aspect ratio as 


^near ” ^0 * A^l^*** 

the lowest order term is found to be a solution of the two-dimensional Laplace equation and 
can be shown to be 

n* - _ p(z* V 1 4* 

^near “ °'^b* ^b' cosh q + cos (jj 

where (h,4*) are2elliptic coordinates centered about mid-chord (see fig. 2). The near field 
solution to 0(A”^) satisfies the equation 


3^P'i 

3n' 


near 


aV 


= -^ g''(z*, y jj) (cosh n sin (() - i sin 2<J)) 
34* A 


the general solutionof this equation is of the form 
P’near '^b> coshtAos T 


+ ^ nq sin n(|) 

n=l 


1 . * . ,1 1 . 

+ g''(Zjj» 'i'lj) (2 sinh q sin (j) + g sin 2<{) ) ( 4 ) 

A 

In the far field region, the characteristic length may be tak^ to be of the order of 
the span in all directions. It is shown in reference 17 that, to 0(A”^), the far field solution 
may be represented by a line of dipoles located along the quarter-chord line of the blade. 

P far “ Pdip^'"’ ^ ’ ^b’ '^b^ 

It is shown in Appendix B that the field of a line of dipoles can be conveniently expressed in 
terms of prolate spheroidal coordinates (see fig. 3) as 
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( 5 ) 


00 

n=l 

The near and far field solutions are then matched by requiring that 

Lim p' = Lim p', 
near far 


L-oo 

b 


r 

— ► o 
s 


which ensures that the behavior of the near field solution at distances of the order of the 
span is like that of the far field solution at distances of the order of the chord. Such a 
matching is achieved by choosing 

h(z^»V = S<V'*'b^ 

where T, 6 x' correspond to r = |, y = ^, = z^. 

V = Avrr;? w pK’ • 

b n=l 

The composite solution for the singular part is obtained as 


p' . = (p* ) . + (p' ) . - (p' ) . 

*^sing near sing far sing common sing 


where 


p' = Lim p* 

^ common near 

E"" 


= Urn p',^ . 


The complete singular part is then 


P'sing “ Pdip^ ^ V * S<2b»V2A r/(Rj- Rq) ~ cosh^J, + 


T1 + COS (j> 


+2(cosh n sin ^ -^sin 
2 

+ — 2 Ijl n sin <|» + sin 2 (j) - ^sin X 


dzj^*- 2A 
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Regular solution. -• The regular part of the pressure field must satisfy the following 
pressure gradient condition. 


+ 




As observed in the previous section, the near field sol^Jtion to lowest order is a solution of 
the two-dimensional Laplace equation which ^ 0(A"^). Consideration of the equation for 
the next higher order only adds terms of 0(A~’^) and can therefore be neglected. It is shown 
(ref. 17) that the near field regular solution is 


^■^2 <P'„ear>rag= f 2 <V P3»bj2s(‘- * 

-1 -2 
The first term, of 0(A" ) 2 can be shown to have a far field expansion of 0(A~ ) and hence 

the second term, of 0(A~ ) need not be considered for the far field. By matching the near 

and far field solutions cuid combining them with a common part, it is shown in reference 17 

that the complete regular part of the pressure field is given by 
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|f j( %) <■ Fj( Vjjjj 2 ^ (*■ Rj) *2r/bj 

+ rp Of F (f >1— (l sin X cosh sin 
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sinh y 
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sin y sin 6 
(cosh M'-cos 0)sinh Y 


(7) 


The complete solution for the pressure field can now be written as p' = p' . + p' , 

and is completely known except for the function g (zg, '?. ) that occurs in the smg§Iar pa^. 
To determine this function, the transverse pressure gradient on a fluid particle is 
integrated along a linearized trajectory from far upstream to a point on the rotor disk. 


9 



This yields the transverse velocity induced by the pressure field at that point, which can 
then be set equal to the known transverse velocity on the blade at the same point. Since the 
unknown function occurs inside the integral, the normal velocity boundary condition leads 
to an integral equation of the form (eq. A^) 




( 8 ) 


The expression for the known velocity on the blade and the form of the linearized 
trajectory are presented in Appendix A. 

Computational Scheme 

The standard method of solving the above integral equation (eq. 8) for the unknown 
function is to assume a collocation function for g (zg, with a finite number of unknown 
coefficients and then satisfy the equation at the same number of collocation points to 
generate a system of equations for the coefficients. Since the unknown depends on both 
radius and azimuth, it would have to be represented by a finite sum of products of radial 
modes and azimuthal harmonics. Van Holten (ref. 17) assumes the following form 


g V = 
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sin kV 



(9) 


This involves a total of (N+1) (2K+1) unknown coefficients. The radial modes are dictated 
by the result of matching and the fact that the dipole line field has been expressed in 
prolate spheroidal coordinates. The field of the dipole line now has the form 

Pdip^'*'’ = “sinh vfcosh V - cos 6) [^^00 E j^Ok *^'^b'*’ ®0k ’^'^bj j 

( 10 ) 

" ^ ko " E i^k *<'^b " ^nk^*" ’^'^bj 1 

n=l L k=i( f-i 


The first term of the collocation function has been used to represent the contribution of 
the regular part of the solution to the far field. 

The pressure field of the dipole line given by equation (10) can be substituted into the 
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boundary condition of equation (8). The resuiting integral equation is then evaluated at each 
of the collocation points which are located at Tu > Tu • The collocation indicies r and s are 
respectively associated with the radial modes af>d alimuthal harmonics. The result is a 
system of linear simultaneous equations for the unknown coefficients that can be written 
symbolically in the form. 


Wfa * * 

^•’br^'^bs^ - '^rs = ^s,oo ^00 


^ ^ ^ rs,okc ^Ok * ^^rs,oks 
k=l 






w 


rs,nkc 


^nk ^ "^rs,nks 



( 11 ) 


Since 1 <r< (N+1) and l<s<(2K+l), there are (N+1)(2K+1) such equations, representing the 
boundary condition applied at that many points on the rotor disk, to solve for the same 
number of coefficients. 

Each of the starred quantities represents an azimuthwise integral, from -<x>to , of 
the function associated with the corresponding A., or B... Since the pressure fielcP also 
contains terms that are completely specified by^W regular part of the solution, the 
integral contribution of these terms appears as w* . The first term of equation (11) is the 
specified normal velocity on the blade due to the*i?lade motion. The starred terms may be 
considered to be influence coefficients since they represent the contributions to the 
induced velocity of those portions of the pressure field corresponding to a unit value of 
their associated coefficients. 

Mention must be made of the behavior of the integrands in the azimuthwise 
integrations. The field of the dipole line becomes singular at its origin, namely, the 
quarter-chord line of the blade. However, this singularity is cancelled out by an equal and 
opposite one in the common part, leaving only the near field solution in the vicinity of the 
blade. Further, the near field solution itseif has a singularity at the ieading edge of the 
biade, of the thin airfoil type (square root). This singularity is integrabie and is evaluated 
by numerical integration up to a point close to the leading edge and replacing the 
remainder of the integral with a simplified analytical expression. Complete expressions for 
the starred integral coefficients are written out in Appendix C, for a singie biade. For the 
case of a’ multi-bladed rotor, the computation is simply repeated successively for the 
pressure fields of the other blades by assuming that the pressure field of each blade to be 
identical in form but shifted appropriately in azimuth dependence relative to the reference 
blade. 

The integrands also require the computation of associated Legendre functions of the 
first and second kind, P j.(x) and Q (x). Both may be computed by using the recursive 
relations for solutions of tiie associated Legendre equations. However, the function of the 
second kind decreases rapidly with an increase in the argument and, beyond a certain range, 
use of the recursive relations leads to an unacceptable loss of significant figures. In this 
range of the argument, the functions are computed using asymptotic expansions. The 
scheme of computation used for these functions and explicit expressions for some of them 
are given in Appendix D. 

The numerical integration utilizes the Gauss-Chebyshev scheme by dividing the entire 
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azimuth interval into sub-intervals of size and applying the rule over each. 


where 
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Although the integration is supposed to begin at an infinite distance upstream, the 
computations are begun at a position where the fluid particle is about three rotor radii 
upstream of the collocation point. In accordance with the form chosen for the collocation 
function, (N+1) points are chosen along the blade span at each of (2K+1) azimuth locations. 

The input required for the computation can be classified under the following headings. 

(1) Blade geometry; root radius, aspect ratio, number of blades, linear twist 

(2) Flight condition: forward speed, inclination of tip path plane to flight path 

(3) Blade motion; collective pitch angle at the root, coning angle, cyclic pitch 
coefficients in the tip path plane (or, equivalently, flapping coefficients in the control 
plane) 

Due to the assumption of linearized potential flow, it may be expected that computed 
airloads would differ from measured airloads because of differences in the lift curve slope. 
Neither can the computed loads be expected to satisfy flapping equilibrium about the hinge, 
as measured loads would. Since the solution features that are of primary interest are the 
variations (chordwise, spanwise and azimuthwise), the use of specified blade motion 
parameters (taken from experimental measurements) would only obscure the comparison of 
predicted variations with measured values. For this reason, it is better to consider the 
blade motion parameters eis unknowns and solve for them using appropriate additional 
equations. If only first harmonic pitch and flapping are retained, then the additional 
unknowns are four in number; collective pitch, coning angle and the two first harmonic 
flapping coefficients. The four extra equations are; 

(1) Average value of total blade lift x number of blades = given rotor thrust 

(2) to W Moment equilibrium of a rigid blade about the flapping hinge (constant, 
first harmonic cosine and sine components). For the flapping hinge on the rotor hub, the 
constant component is determined by the blade inertia (Lock number), while the first 
harmonic components are both zero. 

This procedure is used in all the computations of this investigation. The final 
equations are presented in Appendix E. The augmented set of equations, (N+1) (2K+1) + U 
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in number, are solved simultaneously for the (N+1) (2K+1) unknown coefficients in the 
pressure solution and the four blade motion parameters. With the pressure solution 
completely known, auxiliary quantities such as sectional lift, sectional pitching moment, 
sectional center of pressure, total blade lift, aerodynamic moment about the hub and 
center of the sectional lift distribution can be readily computed. Expressions for these 
quantities are also given in Appendix E. 

Other Methods Used for Comparison 

Wilmer's method . - This method (ref. 6) is an attempt to extend the Prandtl lifting 
line model to rotary wings. The blade is modeled by a lifting line along the quarter-chord. 
The azimuthal variation of the bound circulation is assumed to be quasi-steady, so that the 
wake consists of only trailing vorticity. That part of the wake immediately behind the blade 
is modeled by a plane, semi-infinite trailing vortex sheet so that the velocity induced by 
this sheet is obtained from simple lifting line theory for fixed wings. The rest of the helical 
wake is replaced by a series of infinite, plane vortex sheets under the blade. The induced 
velocity contribution of each layer is given by a spanwise integral that can be evaluated by 
either simple numerical integration or the theorem of residues. The distribution of bound 
circulation can then be determined by a collocation method. The entire computation is 
simple and fast. The solution may be expected to be bad near the downstream edge of the 
disk where the wake curvature is large, thus the assumption of plane wakes is not valid. 

Actuator disk approach . - This model can be formulated by starting with a lifting line 
for the blade and rigid helical sheets for the wake. The induced velocity contribution of the 
wake from a blade at any field point will then be given by an integral, the result of which 
depends on the coordinates of the field point and the aximuth position of the blade (ref. 2). 
The simplification consists of averaging this result over the azimuth, which is equivalent to 
considering the single blade (or finite number of blades) to be spread out over an infinite 
number of blades, each carrying an infinitesimal load. The discrete helical wake is also 
averaged, becoming a skewed, semi-infinite cylinder filled with vorticity. Since the wake 
is a surface of velocity potential discontinuity, a disk section of the wake cylinder is 
replaced by a disk of flow doublets. The induced velocity contribution of this disk is 
integrated over the length of the cylinder to obtain the total wake effect. The 
computational scheme for this model, as described in reference 2, reduces the calculation 
to a series of simple operations, through the use of special functions and integral 
representations. The basic assumptions of the actuator disk approach may be expected to 
be less acceptable at low forward speeds where a discrete, distorted wake representation is 
necessary. 

Segmented vortex wake method . - This approach also starts with a lifting line model 
for the blade and a wake representation as a helical sheet. The simplification consists of 
discretizing the continuous variation of flow quantities on the blade and in the wake. The 
bound circulation distribution on the blade is approximated by a stepped variation which, in 
turn, leads to a trailing wake made up of a finite number of line vortices each carrying a 
finite circulation. The continuous rotation of the blade is also discretized into impulsive 
movement between a finite number of azimuth stations. This leads to a shed wake 
consisting of a finite number of shed line vortices carrying finite shed circulation. The 
result is that the continuous wake sheet is replaced by a vortex lattice of straight line 
segments, streamwise and radial. The induced velocity contribution of each segment is 
known directly as a function of its position relative to the blade. The total wake effect is 
then obtained by a finite summation so that the computational scheme is quickly executed. 
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The disadvantage of such a model, as it stands, is that unsteady effects are not adequately 
represented without a large number of azimuth steps. On the other hand, the approach has 
the advantage that wake distortions can be handled readily so that it may be used either for 
an iterative calculation of a distorted wake or to accommodate a prescribed, distorted 
wake from experimental observation. 

For all the methods considered for comparison, the blade motion parameters are 
considered unknown and solved for in the same way as in the computational scheme for the 
asymptotic method. 


RESULTS OF THE STUDY 
Analytical Study 

In the course of developing a solution for the pressure field of a rotor blade in 
forward flight, Van Holten (ref. 17) considers some simpler flow problems. The first is an 
application of the acceleration potential for the analysis of steady flow past a two- 
dimensional thin airfoil for which the classical solution is obtained. Next the acceleration 
potential formulation, together with the matched asymptotic expansion technique, is 
applied to the problenvof steady flow past a straight, rectangular wing, for which solutions 
to 0(1/A) and 0(1/A^) are obtained. Van Holten shows that the soli^tion to 0(1/A) is 
identical with Prandtl's lifting line solution while the solution to 0(1/A'^) is identical with 
the extended lifting line solution of Weissinger (ref. 18). 

In the current study, the acceleration potential method has been applied to the 
problem of a two-dimensional thin airfoil with oscillatory downwash, as described in 
Appendix F. Particular solutions were obtained for two cases: an airfoil in harmonic 
pitching and heaving motion, and a stationary airfoil with a harmonic vertical gust 
superimposed on the steady stream. These solutions are identical with the results of the 
velocity potential formulation. This is to be expected since the velocity and acceleration 
potential approaches are completely equivalent under the assumption of small disturbances. 
However, the exercise presented in Appendix F does verify the expressions and procedure 
used for the unsteady problem. 

The method was also applied to the problem of a straight, rectangular wing with 
oscillatory downwash, as described in Appendix G. The result of the analysis is an integral 
equation for the spanwise variation of the sectional lift distribution, correct to 0(1/A). An 
equivalent analysis has been presented by Reissner (ref. 19), based on a vorticity 
distribution, the result of which is an integral equation for the bound circulation. 
Numerically calculated lift distributions from the asymptotic method are compared with 
Reissner's results for the harmonically pitching wings presented in reference 20. There is 
good agreement between the two sets of results. 

These two exercises represent an extension of Van Holten's analysis to the case of 
oscillatory downwash, and the results confirm the validity of the formulation and the 
detailed expressions used. It must be mentioned that, in the case of the finite wing, the 
reduced frequency of oscillation is assumed to be of 0(1). Van Holten further applies the 
pressure method to the case of a harmonically pitching wing in yawed flow and points out 
the approximation involved in using a simple cosine sweep correction. 

The formulation of reference 17 for a rotor blade would therefore seem to be 
applicable to the unsteady, three-dimensional^ yawed flow environment of a helicopter 
rotor blade in forward flight, correct to Od/A"^), under the restrictions of incompressible, 
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potential flow and small disturbances in flow velocity relative to forward speed. However, 
the results of such an analysis would be applicable to the rotor blade only to the extent that 
the blade flow environment is itself within the above restrictions. The assumption of small 
disturbances implies that the wake of the blade has a rigid helical shape with all its 
elements being convected with the stream velocity relative to the blade. For this to be 
acceptable, the freestream component of the rotor through-flow must be much greater 
than the rotor-induced component. Under conditions of low inflow, which occur at low 
forward speeds and low thrust, the wake is distorted from the helical shape by interaction 
with the reference blade and the following ones, and also by self-interaction between the 
wake elements. It may be expected that linearized theory will begin to break down when 
these conditions are approached. 

The following factors can also have an influence on the rotor flow field. 

(1) The rapid roll-up of the vortex sheet near the blade tip can form a strong, 
concentrated vortex tube. While the wake of a fixed wing also rolls up at some distance 
behind the wing, the greater concentration of load near the tip of a helicopter blade causes 
the roll-up to be very rapid. In turn, the tip vortex influences the flow near the blade tip 
and also has an influence on the following blade, if it should pass close to it. 

(2) The combination of forward speed and rotational motion can lead to a region of 
reversed flow on the retreating side of the blade. It may be noted that this region grows 
with increasing forward speed. 

(3) Blade pitch angles are varied over the azimuth in such a way that they are 
larger on the retreating side, especially near the root for twisted blades. This can lead to 
conditions of dynamic stall which, like reversed flow, can only be accounted for empirically 
at present. 

(4) The superposition of forward and rotational motion also leads to a component of 
flow along the blade span, which can influence the sectional properties through its effect 
on the boundary layer. This radial flow is strongest along the fore-and-aft diameter of the 
disk and increases with advance ratio. It has been observed that linear thieory computations, 
even under otherwise favorable conditions, are worst at the 0° and 180° azimuth positions 
(ref. 13), which may be at least partly due to this effect. 

(5) At high blade section Mach numbers, compressibility effects and the 
possibilities of local shock formation have to be considered. 

In view of the above discussion, it would appear that any linear theory, such as the 
asymptotic approach, would be applicable only in a restricted range of flight conditions. 
This range of validity can be made more precise only by comparison with appropriate 
experimental measurements. 

Comparison with Experiment and Other Methods 

To evaluate the computational results obtained from the asymptotic method, 
conditions from the following experimental investigations were analyzed. (See Table I for 
specific geometric and flight conditions.) 

Case 1; A two-bladed, teetering, model rotor of aspect ratio 5.^ with untwisted, constant 
chord blades was tested in a wind tunnel at advance ratios, y, of 0.08, 0.15 and 0.29. 
Results for this experiment are presented in reference 21. 

Case 2; A four-bladed, articulated, full-scale rotor of aspect ratio 17.2 with twisted blades 
was tested in flight at advance ratios of 0.06, 0.13 and 0.29. The experimental results are 
tabulated in reference 22. 

Case 3; The four-bladed rotor of Case 2 was cilso tested in a wind tunnel at advance ratios 
of 0.29, 0.39 and 0.45. These results are presented in reference 23. 


15 



TABLE I. - GEOMETRIC AND FLIGHT CONDITIONS FOR THE 
EXPERIMENTAL CASES 
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These cases were also analyzed using Willmer's method, the actuator disk method and 
the segmented vortex wake (SVW) method. In all of the computations, it was assumed that a 
fluid element is convected normal to the disk with the sum of the freestream inflow 
component and a constant thrust-induced velocity, determined from simple momentum 
relations. It was also assumed that the rigid blade motion parameters (collective pitch, 
coning angle and cyclic pitch) are unknown. The results presented here include the 
azimuthwise variation of total blade lift, spanwise variation of sectional lift and chordwise 
variation of surface pressure differential. 

Results for the two-bladed, teetering, model rotor (Case 1) are presented in figures 5 
to 11. The total blade lift (see fig. 5) for advance ratios of 0.08 and 0.15 from the 
asymptotic and SVW methods are quite similar and agree fairly well with the measured 
data. For the advance ratio of 0.29 the two methods yield different results and neither 
compares well with the experimental lift. The actuator disk method yields fair agreement 
with the measured data at the higher advance ratios of 0.15 and 0.29 but shows marked 
deviations at 0.08. Willmer's method shows the poorest correlation with the measured 
values. This disagreement appears to increase as the advance ratio is decreased. In 
summary, the asymptotic method provides the most consistent results for Case 1 except for 
the advancing blade (0 < < 180 ) at an advance ratio of 0.29. This deviation has been 

reduced by decreasing the number of azimuthal harmonics and will be discussed later. 

The spanwise distributions of sectional lift are shown in figures 6 to 8 for all three 
advance ratios at intervals of 30° in azimuth. For clarity, only the results of the 
asymptotic method are shown for comparison. Whenever the experimental curves are 
smooth, the computed values have nearly the same shape, except near the tip. At some 
azimuth locations, the measured curves exhibit sharp variations that are not illustrated by 
the computed load. These kinks are probably due to a close encounter with the tip vortex 
from the preceding blade, especially since they generally occur near the extreme advancing 
and retreating blade positions and are most predominant at the lowest advance ratio of 
0.08. 

The chordwise variation of surface pressure differential is shown in figures 9 to 1 1 for 
the advance ratio of 0.29 at radial stations, r/R. = 0.75,0.85,0.95, and azimuthal intervals 
of . The measured variations are seen to be^fairly smooth throughout, with the same 
general shape as a two-dimensional distribution of the asymptotic method. Computed 
results are not corrected for the leading edge singularity but agree with the measured data. 

Flight-test results for the four-bladed, articulated, full-scale rotor (Case 2) are 
presented in figures 12 to 18. None of the Case 2 results from Willmer's method are 
presented, because they differed from the measured data by nearly an order of magnitude. 
The total blade lift for the advance ratio of 0.06 (see fig. 12(a)) as predicted by the other 
three methods shows significant disagreement with the experimental results. This can be 
attributed to the low advance ratio for which the actuator disk method loses its validity 
and the small value of wake spacing. The wake spacing as listed in Table I was computed 
from the total inflow velocity, w^, which is the sum of the freestream inflow and the 
thrust-induced inflow. At this condition of small wake spacing it can be anticipated that 
the effects of wake distortion have caused the asymptotic and SVW methods to yield poor 
results. These effects could also be the reason for the complete failure of Willmer's 
method. As the advance ratio is increased to 0.13 (see fig. 12(b)) both the disk and SVW 
metjiods yie|d bette^ correlation except for the downstream portion of the rotor disk from 
300 to 360 and 0 to 90 . To obtain the asymptotic results for this case of very small 
wake spacing it was necessary to reduce the integration interval from AY = 0.3 (as used for 
the preceding test conditions) to AT = 0.1. Even with this refinement it is apparent that the 
correlation with the measured data is unacceptable. It can therefore be concluded that the 
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asymptotic method yields unreliable results when the successive wake spacing is less than 
five percent of the blade radius. At the highest advance ratio of 0.29 (see fig. 12(c)) all 
three methods provide consistent results for the retreating blade portion of the rotor disk, 
and their correlation with experiment may be considered fair. For the advancing blade all 
methods yield different results and none is close to the measured lift. As previously 
mentioned for Case 1, the asymptotic results for the advancing blade can be improved by 
decreasing the number of azimuthal harmonics and will be discussed later. 

The results for the spanwise lift distributions which are presented in figures 13 to 15 
show the same trends as for the two-bladed configuration of Case 1. The experimental 
curves show sharp variations in the neighborhood of 90° and 270° azimuth. As before, these 
are probably due to tip vortex encounters, which have been discussed in detail by Scheiman 
and Ludi in reference 2k. 

The measured chordwise distributions of figures 16 to 18 for the advance ratio of 0.29 
have an approximately normal appearance at r/R. = 0.75 and 0.85, and correlation with the 
computed values is about the same as for the two-bladed rotor of Case 1. At r/R, = 0.95, 
however, the mecisured data show abrupt variations, again around the 90° azimuth ^sition, 
which are not apparent in the computed points. It may be mentioned here that measured 
chordwise distributions have been compared with static two-dimensional data (from wind 
tunnel tests) for this rotor in reference 25. That comparison indicates that the measured 
data have a two-dimensional appearance over most of the disk. 

Wind-tunnel results for the four-bladed, articulated, full-scale rotor (Case 3) are 
presented in figure 19. As for Case 2 none of the results from Willmer's method are 
included due to the poor correlation with test data. Both the actuator disk and SVW 
methods yield comparable results and agree fairly well with the experimental data at all 
three advance ratios except for the downstream portion of the rotor disk 020° to 360° and 
0 to 90 ) at p = 0.45. All asymptotic results are based on the reduced integration interval 
of 0.1. This method shows acceptable agreement for the advance ratios of 0.29 and 0.45 
except again for the downstream portion of the disk at p = 0.45. For the advance ratio of 
0.39 the asymptotic method yields unacceptable results. It will be seen later that the 
correlation for this condition is significantly improved by decreasing the number of 
azimuthal harmonics. 

In summary it may be noted that all methods of analysis provide less than desirable 
comparisons for all three Cases at the highest advance ratios for the downstream portion of 
the rotor disk. It is likely that at high advance ratios the downstream portion of the disk 
can be affected by blade root vorticies and the wake of the hub and other fittings. Since 
none of these methods include these effects, such discrepancies should be anticipated. It 
can also be noted that the asymptotic method acknowledges a continuous wake, while the 
SVW and actuator disk methods correspond to lumped and averaged wake representations 
respectively. For this reason the sensitivity of each method to wake distortions may be 
expected to be different. 

It is seen that Willmer's method is most sensitive to low inflow conditions. For wake 
spacing less than approximately 0.06 the results are unacceptable. The segmented vortex 
wake method seems to be less sensitive to these conditions than the asymptotic method and 
this feature could be a consequence of the discretized nature of the method. The actuator 
disk method is not sensitive to inflow conditions since instantaneous wake effects are 
averaged out; however, the results indicate a limit of acceptability in terms of advance 
ratio, below which the total lift is unsatisfactory. This lower limit appears to be 
approximately 0.1. Ormiston (ref. 27) observed a similar limit for the actuator disk method 
to be around 0.15. 
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In the context of these comparisons, it should be noted that a similar study was made 
by Hille in reference 26 for Case 2 at an advance ratio 0.13. The blade was represented by 
a lifting line along the quarter chord, and the boundary condition was satisfied along the 
three-quarter chord line. The method used a continuous vortex wake trailed from the 
instantaneous flapped position of the blade. In spite of refinements such as using a linear 
variation across the disk for the thrust-induced component of the total inflow velocity, it is 
noted that the comparison is not satisfactory. 


Use of Measured Section Characteristics 

The possibility of using measured airfoil section properties in the asymptotic method 
is discussed by Van Holten in reference 28. Since the near field asymptotic solution, to 
lowest order, is that for the pressure field of a two-dimensional airfoil, it would seem 
possible to introduce measured characteristics in this part of the solution, to achieve better 
correlation with measurements. However, a part of the near field solution, the function 
g(z. ,'i'j^) is undetermined by itself. It is determined only by combining the near field, far 
field Sid common part to form a composite solution. It is then made determinate by 
applying an essentially three-dimensional boundary condition, namely, integrating the 
pressure gradient on a fluid particle as it travels along a helical trajectory. For a complete 
two-dimensional correspondence, the path relative to the blade should be a straight line 
from an infinite distance upstream, as in flow peist a fixed wing. 

Section characteristics are measured on a two-dimensional airfoil in a wind tunnel, 
with the airfoil being given a pitching and/or heaving transverse motion. The form of the 
boundary condition for a rotor blade section (Appendix B) shows that not all of the terms 
can be simulated in this way. The direct use of wind tunnel data in the form of empirical 
factors in the asymptotic approach requires careful consideration of these boundary 
conditions. 

Computational Study 

As previously discussed, the unknown function in the pressure solution, g (zu>T. )> is 
determined by using a collocation technique which includes a sum of functions, each &eing 
the product of a radial mode and an azimuthal harmonic, multiplied by a constant 
coefficient. With this representation, the solution requires the determination of these 
coefficients by forming a system of linear equations to be solved simultaneously. If five 
radial modes and five azimuthal harmonics are used, as was the case for all the 
computations previously discussed, there are a total of 55 unknown coefficients and the 
system of equations has a matrix of size (55 x 55). Since the determination of each 
coefficient matrix element involves an azimuthwise numerical integration, the total 
computational time is determined primarily by the time required to set up each coefficient 
element and the number of such elements involved. The collocation form determines the 
latter while the efficiency of numerical integration decides the former. To study the 
efficiency of the computational scheme, it weis decided to vary both of these factors. 

The numerical integration uses a five-point Gaussian formula over an azimuthal sub- 
interval, AT. For the previous computations this interval size was chosen to be 0.3 (rad) for 
all of Case 1 and for advance ratios of 0.06 and 0.29 for Case 2. The other conditions all 
required AT = 0.1. The computation time increases threefold with a reduction in AT by the 
same factor. Other conditions studied with AT = 0.1 were Cases 1 and 2 at an advance ratio 
of 0.29. The results are illustrated in figures 20 and 21 from which it can be seen that the 
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refinement makes no significant difference in the total lift at these conditions. In an 
effort to improve the poor correlation of the asymptotic method for Case 3 at the advance 
ratio of 0.39, the integration interval was reduced from 0.1 to 0.05. The resulting lift as a 
function of azimuth angle showed no significant improvement as illustrated in figure 22. 

The number of radial modes was increased from 5 to 10. This doubles the number of 
unknown coefficients and approximately quadruples the computation time. This 
modification was also only applied to Cases 1 and 2 with y = 0.29. For Case 1 the results 
were only marginally different. However, for Case 2, the results were unacceptable which 
may be an indication that, when the method is applied near its limit of validity, an increase 
in the number of radial modes may make the computational scheme unstable. Reducing the 
number of modes to four led to unsatisfactory results in some cases where the five mode 
computation had yielded acceptable results. It would therefore seem advisable to set the 
number of radial modes at five. 

Varying the number of azimuthal harmonics does not seem to affect the solution as 
much as the number of radial modes. Computation with five modes and three harmonics for 
the conditions of Cases 1 and 2 showed that in most cases there was no significant 
difference from the five and five computation. In fact, for the advance ratio of 0.29 in 
Cases 1 and 2, the five and three computation led to significantly better results, as can be 
seen in figures 20 and 21. The correlation with the measured data is greatly improved for 
the advancing blade with the five and three computation. The reason for this improvement 
is not clear. It is possible that a method based on linear theory is not well suited to the 
calculation of higher harmonic variations so that it is better to neglect them altogether 
than to include them and miscalculate their contributions. The five mode and three 
harmonic representation was also used for Case 3 at the advance ratio of 0.39. This 
resulted in a most dramatic improvement as illustrated in figure 22. Based on this limited 
study, a computational scheme based on five radial modes and three azimuthal harmonics 
would appear to be an optimum compromise between accuracy and computational expense. 

The computer program has been checked out and found to be free from programming 
errors although it could possibly be made more efficient with respect to computer time. As 
it stands, for the two-bladed rotor at an advance ratio of 0.29, with the five and three 
scheme, the program requires about 1.5 minutes to execute on the CDC Cyber 70 Model 74- 
28 computer. For different conditions, this time would scale approximately in direct 
proportion to the number of blades and in inverse proportion to the advance ratio. As a 
comparison of computation times, it may be noted that the asymptotic computation scheme 
with five radial modes and three azimuthal harmonics required approximately six times as 
much execution time as the segmented vortex wake scheme with five spanwise and twelve 
azimuthwise segments. The accuracy of the results was comparable for these two schemes. 
However, the computation time depends strongly on the size of the system of equations 
involved in the method. As an example, the computation times for these two schemes would 
be about the same if the vortex wake method used nine spanwise and sixteen azimuthwise 
segments. 


CONCLUDING REMARKS 

The combination of the acceleration potential formulation with the matched 
asymptotic expansion technique leads to a systematic determination of airloads on a 
helicopter rotor in forward flight. For a straight wing in steady flow, the method reduces to 
the simple lifting line solution of Prandtl or the extended lifting line solution of Weissinger, 
depending on the asymptotic order of terms retained. Numerical results for a straight wing 
in oscillatory motion compare well with results from Reissner's unsteady lifting line theory. 
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Under the assumptions of incompressible, potential flow with small disturbances, the 
asymptotic method appears valid for application to a helicopter rotor blade in forward 
flight. However, the rotor flow field is significantly influenced by nonlinear and real fluid 
effects such as wake distortion, vortex sheet roll-up, reverse flow, dynamic stall, radial 
flow, finite vortex core radius etc. In any linear, potential flow analysis, these factors are 
beyond the scope of the basic assumptions and may be expected to cause deviations 
between computed results and measurements. 

The asymptotic scheme was used to compute airloads for comparison with 
experimental results corresponding to a two-bladed, teetering model rotor of medium 
aspect ratio and a f our-bladed, articulated full-scale rotor of large aspect ratio. In general, 
the computed results were acceptable for those cases involving a successive wake spacing 
greater than five percent of the blade radius. Discrepancies at closer wake spacings (lower 
inflow) are to be expected from any linear, potential flow analysis. 

A collocation scheme for the asymptotic method of five radial modes and three 
azimuthcil harmonics appears to be an optimum compromise between accuracy and 
computational expense. 

For further comparison, the same cases were also computed using other approximate 
methods. Willmer's method was unacceptable for all cases involving the f our-bladed rotor. 
The actuator disk approximation showed no sensitivity to inflow and led to fairly 
acceptable results except at advance ratios below 0.1. A method based on a lifting line with 
a rigid, segmented vortex wake compared well with the asymptotic method and showed 
lesser sensitivity to low inflow conditions. 

The near field solution for the rotor blade, to lowest order, involves the solution of a 
two-dimensional Laplace equation and therefore has the capability of using measured 
airfoil characteristics. However, the form of the boundary condition on the blade is 
essentially three-dimensional, and the use of measurements made on two-dimensional wind 
tunnel models requires careful consideration of these boundary conditions. 

A comparison of computational times between the asymptotic method and the 
segmented vortex wake model indicated that the asymptotic approach required 
approximately six times as much execution time as the vortex model. This comparison was 
for comparable accuracies. It should be noted that for both techniques the execution time 
is highly dependent on the number of radial and azimuthal collocation parameters. 
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APPENDIX A 

ROTOR BLADE BOUNDARY CONDITION 


Each blade of the rotor is presumed to be rigid, uncambered and rectangular in 
planform with linear twist. Coordinate systems which are used to describe the rotor 
configuration are illustrated in figure 1. The "rotor-axes" coordinates, x , y , z , are 
oriented such that the x -y plane is the tip-path plane and the freestream direction is 
parallel to the x — z plane- *The "flow-axes", x, y, z, are inclined to the rotor-axes by the 
angle a , which is a rotation about y , such that the freestream velocity is in the negative x 
direction. A third system which is fixed to the blade has its origin at the quarter-chord of 
the mid-span. These "blade-axes", Xj^, yj^, z^, are inclined to the tip-path plane by the 
coning angle, a , such that coincides witn the quarter-chord line. The rotor angular 
velocity is assumed constant. In the blade axes system a general point on the blade is given 
by 


yb= {Qr-^^'b- ""b 


where 0 is the pitch angle at the root, relative to the x. - z. plane. The total twist, e , is 
positive ^or section incidence decreasing toward the blade tip. The flapping angle is 

3 = a^ - aj cos T^j-b^sin 

so that the root pitch angle is given by 

0 = 0 + b , cos 'P. - a , sin Y. 
r 0.1 b 1 b 

Consistent with linearized theory, small angle assumptions are made so that a point on the 
blade surface is given by 

X = -rj^ cos Y|^ - Xjj sin 
y = -xj^cos Y^ + r^sin Y^ 

2 = {0Q ■ ^^"b - - M ""b - ^"b ■ ^o 

+ X + (bj Xj^) cos Yj^ - (a^ x^^) sin Y|^ 

The velocity components at this point are given by 

u = ^ = -( ft Xj^ + Tjj) cos Yjj + ( ft rjj - Xj^) sin Y|j 
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V = ^ = (fl fjj - Xjj) cos 'i'l^ + ( + r^) sin 

w = §f =-«bjX^sin 'Fj^-flajXjjCos 'i'b + x^{6o- e( - Ro>/(Ri - Rq) 

+ bjCos Yjj-ajsin 'i' 5 } + {-e Xj^/(R^ - R^) - a^} -Ua^ 

• • 

Assuming u = U + u', v = v', w = w', the quantities Xj^ and Tj^ can be determined as 

Xfa = fl r + U sin - v' cos 'F - u' sin y 

fjj = U cos Yjj - ftXjj - u' cos ^b + V sin 


Substituting these values into the expression for w, linearizing and setting y = U/j2R,, the 
normal velocity of a point on the quarter-chord line is obtained as 


(wjj/flRj) _ pe^ “ s ^^’b ■ ‘‘b^^i ■ v“r “ 

b 

+ [bi r^/Rj - ya^] cos + [- aj r^/Rj 

+ u{eo- e(rb-Ro)/(Ri“ VH '*'b 

+ [y aj/2j cos 2 + [ybj/2] sin 2 (Al) 


It has been assumed in the foregoing that the rotor incidence, a , and the coning angle, a , 
are small. *" ° 

It can similarly be shown that the normal acceleration of a fluid particle on the blade 
is given by 




Dw^ 


£)-j ~ ■’k/Ri) x,^/R| + F,( Yk) + F,(Y ^,) r^^/R 




2' "b' 




(A2) 


where 
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e Rq/(Rj - Rq) + 4e (r^/Rj) [RjARj - Rq)] 


■ ®o ■ 


- 2bj cos 'Fj^ + [2a j + 4 ye Rj/(Rj - Rq)J sin 

^2^ ' \} ^ ^ U eRfl/^Rl “ RqO \ ^ \ 

+ 2 ybj cos 2 'F|^ + 1^- 2 ya^ - y^e Rj/(Rj - Rq)J sin 2 Yj 
F-^i 'V^) = 3q + jj-2aj - ^ y e Rj/(R^ - Rq)] cos - 2bj sin 'Fj^ 
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With gy- , the solution for the pressure field must be obtained such that the normal 

pressure gracueht on the blade surface balances the normal fluid acceleration. This yields 

Further, the normal velocity boundary condition at the quarter chord must be satisfied by 
integrating the normal pressure gradient along the linearized trajectory up to the reference 
point. A point (x^, y^, z^) on the quarter-chord of the blade at time given by 

’'bo 


V‘^1 = '^bo 


V'l = ^o<‘-''bo'“l>-“rW‘'l>’=“ ’’bo 


The linearized trajectory in the flow axes system is given by 

x(t) = x^-U(t-t ) 
o o 

y(t) = y„ 

z(t) = 

2H 



and, in the blade axes system, by 


WMCfb- V5o)sln fb 

yb/R, = y(Cr+aoCos ^b^'I'b- fbo' ' '•'b ' 'i’bo’-'l 
V^l = - <^0 * Rl> * M< \ - fbo) cos ?b * <Sbo/Rl> < '^b - ^0* 

The normal velocity boundary condition at the quarter-chord of the blade is then written as 

'J'k 

bo 

4<^b„.V = -/ (A*) 


APPENDIX B 


LIFTING LINE OF PRESSURE DIPOLES 

The pressure field must satisfy the Laplace equation which can be written in 
cylindricai coordinates (r, x , z) as 

9^^ iiE > + 3^ ,0 (Bl) 

r 3r >“ 9x 

To apply a separation of variables, let 

p = R(r) X( x) Z(z) 

This leads to the ordinary differential equations 

Z" + q^Z = 0 
X” + n^X = 0 

r^R" + rR' - (n^ + q^r^)R = 0 


The first two equations have sines and cosines as elementary solutions, while the third has 
modified Bessel functions of the first and second kind, I (qr) and K (qr), as solutions. 

For the solution to represent a lifting line, it must tend to zero as r-*-<» and become 
singular as r-^ 0, which excludes the use of I (qr). A field point can be represented by a given 
value of X or with any multiple of 2 ii added to it, which requires that the solution be 
periodic in X and that n be an integer. Since the soiution must be antisymmetric about x =0 
or IT, only sine solutions can be allowed for X(x). Further, for the lifting line to be built up of 

dipoles, n=l, since K (x)~x~'^, x-^0, n=l, 2... ^ t \ 

The general solution for a dipole line is then of the form, p(r,x»z) = sin x 
{ A(q) cos qz + B(q) sin qt } . To obtain a line of strength f ^(z), - s < z < x, the functions A and 
B are chosen as 


A(q) = q 


1 


B(q) = q 


s 


f ) cos q^ d 5 


f j( sin q ^ d 5 
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Since all positive values of q are allowed, the dipole line has the solution 

Pdip^'’’ X »z) = I t fi(?)dc qKj(qr)cos {q(c -z)} dq 


(B2) 


Integrating by parts and making use of f j(+s) = 0, this expression can be alternately written 


as 


Pdip^‘'’X»2) = -^^ 


JL / (c-z)dc 


1/2 


It can be seen that a field for the lifting line of this form involves numerical integration 
along the span. An alternate approach is to obtain the field in prolate spheroidal coordinates 
which are related to cylindrical coordinates through (see fig. 3) 


r = s sinh Y sin 6 
z = s cosh Y cos 0 


X = X 


The Laplace equation becomes 


sinh" 


i + coth T ^ P- + cot 0 ^ 

y e '' 8 '^ sy 3,2 aej 


sJe.o 

sin^e 3x^ 


(B3) 


Separation of variables as 


p = \p(y)0(0)X(x) 


leads to 

+ (coth T) \1^'- |n(n+l) + — ^ — | \1^ = 0 
1 sinh"^ T ' 

0" + (cot 0)0'+ 1 n(n+l) - 0=0 

' sin^ 0 ' 


X" + m^ X = 0 
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The solution for X(x) is as before^^e first two equations have, as elementary solutions, the 
associated Legendre functions, (cosh T), (coshY ), (cos0 ), (cos0 ). For the 
solution to represent a dipole line it is necessary that m=l. I^r it to vani^ as r-*-® {^-t-oo ) 
and be singular along the lifting line (Y=0), only.QMcosh '!') can be considered for \|;. To 
avoid any other singularities in the field, only Pj!,(cos 0 ) can be considered for 0 . The 
solution for a lifting line of dipoles is then of the form 

Pdip('i'»e»x) = ^^ f A^P^(cos0)Q^(cosh'F) (Bif) 

n=l 

It can be seen that such an expression for the field of a lifting line, conveniently truncated, 
involves no numerical integration edong the span. 
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APPENDIX C 


INTEGRAL COEFFICIENTS FOR INDUCED VELOCITY 

In order to simplify the final form of the integral terms in equation (11), the following 
derivatives will first be defined. 


D _ a r sin Y coshY sin 6 ~ 1 
^ \ Lsinh Y (sinh^ Y + sin^ 0 )J 

_ ^ ^1 r sin^ Y I cosh Y 

Unh^ Y + sin^ 0 ( sinh^ Y + si 


3 2 

cosh Y sin 6 


2 cosh Y sin^ 0 (cosh^ Y +cos^0)i 


(sinh^ Y + sin^ 0 


sin^ 0 sinh^ Y (sinh^ Y + sin^ 0 ) 


2 

cos ^ cosh 


sinh^Y (sinh^ 


:osh Y 

2 

Y + sin 0 ), 


D2 = 




r sin Y sin 0 ~ | 

[sinh Y (cosh Y - cos 0)J 


2R 


R, r • 2 
1 sin Y 

rV [cosh^ - 


/ 2 2 2 2 
sinh Y cos 0 - cosh Y sin 0 

2 2 

cos 0 V sinh Y (cosh Y - cos 0 ) 


. 2 
sin 


1^0 (cosh Y + cos fl ) ^ COS^y 1 

(cosh Y-cos 0)^ ) sinh^Y (cosh Y-cos0)| 


r\ 3 r sin 

3=^(^)[P7^ 




O'J 


cos2y 

(r/Rj)^ 
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D,= 


d 

sin (j) 1 

(yh\ 

Vr,/ 

Lcosh n + cos ({) J 


2AR 
= (R 


:aKj - 
.(si 


.(sinh^ n 


sinh n 

2 ^ 
cos (j) 


cos d) + sinh n cosh r\ cos 2 6 


+ cosh^ n 


sin"^ 6) (cosh n +coS(f))‘' 


-1 


D j ® ^ Qj^fcosh Y )] 

ir) 

R 

= 75— Vt r . Q u P^(cose)Qi(cosh'r) 

(R^-Rq) Lit sin 0 sinh 'F n ^n 

+ — I sin 9 cosh Y pj^(cos e)Q^(coshY ) 

2Ti(sinh^Y +sm 0) ' " " 

- cos 0 sinh Y P^(cos 0 ) Q^(cosh Y)| J 


^6 = 


Hi) 


[■ 


e"^ sin (j)J 


2ARj 

(R^RJ* 


2 2 
e~ ^ (sinh n cos d) - cosh n sin (1>) 

2 2 2 2 
sinh n cos (j) + cosh n sin (j) 


‘] 


It must be noted that the complete expression for the pressure field is valid only for a field 
point with a spanwise coordinate within the blade span. For a field point outside the blade 
span, it is assumed that only the far field exists. In the expressions that follow, that part of 
the pressure field that is valid only within the blade span is multiplied by the factor 
(H(zg+1)-H(zg-1)} . The Heaviside function, H(x), that is unity for x>0 and zero for x<0 will 
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be denoted by H, so that 


H = 


fl, |Zbl<l 

0, |zjj I > 1 


In terms of the derivatives (Dj to Dg), the integral coefficients are as below 


w 


rs,oo 






/ 


+ I HA (1 + 2 jj)D^dYjj 


Ad ♦ 


■'bs) 

V g” '‘'bs 
wsin 


coth( n j/2) 


(^coth(rii/2)- nj) 


where 


cosh 




"bi'^ * >■ ’^bs) *1 
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'^bs 

<s,okc = f [^2' *^ 3 ] •< '^b d 'J'b 




J' HA(1 + z^)D^coskYjjdYjj 


Ad.z,^) 


(^^ysin^'bs) 


cos k 'Fjjg coth( n j/ 2 ) 


pcos k f cos f - 1 ( 1 - ^)(l * 4 ) sin k¥^' 




i^* ysin 


(icoth^-nj) 


k (*- 4“)>- -=« 'f bs '^'^bs^8^(‘-r°n‘^v) 


bs 


+< 


(•^tysinl'J 


/l 1 \ 

I ^ coth ~2 - sinh n j ) 


cos k ¥ 


t ysln 


bs r sin Y*sln 9 * oa/i . * \\ 

\ Lsinh ^ • (cosh Y* - cos 6 ') " br^J 


b IT * * 

where ( Y 0 ', x ') corresponds to r = 2 > x = Zjj = Zjjj. 
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* * 
"rs^ks *rs,okc 

-cosk ^bs 


except replace cos k Y with sin k 'i' and sin kY with 




pX>D 4« 


b 


-X )77 


coth( n j/2) 


(•^ * \s) 


, 0 , » 


PCosTbs "<"*‘>Pn<V 
TT ^ < 2“* T 

■‘“"V 


2ir 


- n,) 


p^cos^'F. n(n+l) P^(zj^*) , Pi 

(icoth^-sinh np 


(^^PSin^bs/ 


J. 

* >* \s) 
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* 

rs,nkc 


V 

/ ■ [°3 ’’nK> >^3] <=“ I' 'fb “ 1 


bs b 


/U 

coth(-^)cosk'l'b3 


-^^f^ p‘(^J 


"br n'^-br'/r 


Msinfb,) 



l< V‘‘V^ ‘’Xr> bs- ‘’Xr> ’’bs '< '*' 

(l"th^-b,)/(^t usin 


“'ll W coskt^^nWDP^Cz^^) 


UirA 


Jl-z*^ 

\ br 




bs 


'> 

*tIsa (‘-r°)^3/^pi<C>'<^=“'‘''bs 

({coth'^ -slnh- n,y(-|^ + U sin 'll Jj 


COS k ^ 


bs 


(-^ t usin 


[k «'> '^ifcosh f ■) * pXr>] 



w* same as except replace cos k with sin k and sin k 

-cosk'i'j,3 


rs 


/[ 

^ 00 L 




(-Jil.. ..(-^5 


F,Do + H-^ ^ (Fo+^ FJD, 


8A 


2 * 3 ""2 


8 a 2 "2 ^Ri 


d'F. 


OS /* 0 


f ° ^F)D 


6‘'*b 



“ 1 

'‘•hr \ / 

^ '^bsj 

(-S1 


SA^ 

(-^ » Fsin 


/3 ‘^11 

•(^-m-e +2® ) 

,/, ^o\ 

2 


F /w iL\ 

32 A^ 

^1 \^bs’ rJ 


with 
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APPENDIX D 


COMPUTATION OF ASSOCIATED LEGENDRE FUNCTIONS 


These functions are independent solutions of the associated Legendre equation and 
obey the following recursive relations (ref. 29). 


"’ll W = idsT - 1) C2<>‘il 


P*^(x) = 2(m-l) 

" Vi^ 


P*JJ'^(x) - (n - m + 2) (n + m - 1) Pjj''^(x) 


and likewise for (x). The functions P*^ (x) can be computed without difficulty using the 
recursive relations ,'^given the explicit expressions for the first few functions. However, the 
functions (x) decrease rapidly with x for x>l and, beyond a certain range of x, use of the 
recursive rePations leads to excessive loss of significant figures due to roundoff error. In 
such a case, these functions must be computed either by using asjsjmptotic expansions in 
inverse powers of the argument or by using the definition of (x) in terms of the 
hypergeometric function. 

In the computations carried out here, it is necessary to compute many of these 
functions very frequently, which requires that each computation be done as quickly as 
possible. For this reason, explicit expressions were used for the range of n and m required, as 
indicated below 


Functions 



P® = (3x^-l)/2 
P° = (5x^-3x)/2 
pj = (35x'^-30x^+3)/8 


36 



m = 1 


p{ = 

P^ = 3x‘^7 

P^ = (1 5x^-3) VI7 12 
pJ = (35x^-15x)Vu?’ /2 

m = 2: 

Pj = 0 

P^ = 3(l-x^) 

P^ = 15x(l-x^) 

P^ = (105x^-15) {l-xh/2 

m = 3: 

Pj = P 2 = 0 

P^ = 15(l-x^)^^^ 

P^ = 105x (1-x^)^^^ 
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Functions Q'^(x), 1<x<3; 


m = 1: 




pi=V7:r[<i^ -5x1 

3 2(x^-l) -1 

n> In"** O5x'*-30x^ t 3) ,„,x^ Ml 

Ql> - ’=< -1 L » *" x-1 g(^2_,, 8 * 24j 


m = 2: 



Q 


2 

2 



+ 


x(3x^-l) ~j 


q 2 = (x^-l) fi|S - 5 ] 

^ L 2 X-1 (^2_j)2 J 


p2 = (,2.„[I105x^ 



(33x^-15x) 

(x^-l) 


+ 


x(35x^-30x^+3) 

^Kx^-D^ 


- 105 
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xl•=^ 



Functions ^ >3 

m = 1: 

Q‘^-(x2.1)gx-^ix-^fx-^ifX-‘»] 

qJ-(x^..)[^x-^^x-^^x-‘“] 


m = 2: 




Q2 = (x^-U 

;|x-^fx-^fx-’] 

qI~<x^-i) 

[fx-^|x-*.|fx-‘»] 

Q^ = (x^-1) 

[Mx-^ix-’] 
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For some of the computational cases, when the fun<^ions were required for greater order 
and degree, the recursive relations were used with being defined in terms of the 

hypergeometric function teyond a certain value of x. (x) was computed recursively 
starting with Pq, Pq, Pj, P^. The functions Q (x) were computed as 


X >3: 



X <3: 



(n-l)/2 


(2n-l) „m 
n ^n-1 



k=l 


(2n-l-»k) 
(l+2k) (n-l^ 


P 


m 

n-l-2k 


(x) 


1 , x-fl 

*2*" FT 


p>). 




n 


m(m-l) (m~2) (3x + 1) 


m(m-l) (m-2) (m-3) x (x^ 


pnr>-*fe) 

n 


This expression is valid for ms^. 
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APPENDIX E 


OUTPUT PARAMETERS AND EQUILIBRIUM EQUATIONS 


The perturbation pressure on the surface of the blade can be obtained by considering 
only the near field solution and setting ti'=0. 


2 d 2 Pb ^ '^b^ 1 +"005 

pft Rj 


<l> 


2 sin (j> 

^ ^ *2S (*■ Rj) F 3 (f|j)j sm <|. 


R 2 

*'l('‘'b’i^) [sin (J) + 1 sin 2 <|)] 


(El) 


By integration with respect to the chord, sectional quantities such as lift and pitching 
moment can be derived. 

Section Lift: 


*. * V Q 








4A 




(E2) 
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Pitching moment about a point x = a (positive nose-down): 


m*(z;'J'b 


) = 


m 


2 k 
pft R, 





(E3) 


Position of center of pressure from the leading edge: 



(E4) 


The sectional lift can be integrated along the span to give the total blade lift and the 
aerodynamic moment about the hub. 


Let 




1 



P^(cos 0') Q^(cosh T') 
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Total blade lift: 






„l^cosk1'b*Bo|<sink'rj,) 




K 


*10 * <*ik “* '' '•'b * ®ik '< V 


■] 




K 


*n0 * <*rt< b * ®nk b> 


* -2 (>- "2< »b> * ^2 (•- >=3< (' * r“) 


'] 


16A' 


(-5;)’[l-.-(-t:)/(-r;)l 


- 2 bj cos Tb+|2aj+4ue^— ^ I sin j 

V^J 


(E5) 
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Moment of lift about the hub: 



(E6) 





Radial position of lift center: 



(E7) 


The additional equations needed to solve for the blade motion parameters are 


1 ) 


J_ 

2ir 



L*(V<f1'b 


TiCy 

“B 


where B is the nutpber of blades. 
Substituting for L 


[-a(‘- 

V Rj/ J'^00^3 it 



*_L 

n=i 

- » 


27tA 

rJ 0 


1 

8A^ 

(‘■Rl) (**R])®o = '=T - 


(E8) 


ii5 



2 ) 



M*( V d 


H?a 
o_ 


Substituting for M 



!(-?■.)(- 1 ;)« 


10 10 



20 2A 




nO 






3) 


2ir 

i / M*(Vcos = 0 

0 


Substituting for M 






(ElO) 
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4 ) 


I M*( 'i'jj) sin = 0 


Substituting for M 





(Ell) 
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APPENDIX F 


AIRFOIL WITH OSCILLATORY DOWNWASH 

The coordinate system for the airfoil is shown in figure 2. Under the assumption of 
small disturbances, the perturbation pressure p satisfies the Laplace equation 


2 2 
2 + 2 - ^ 


Since oscillatory downwash is considered, the complex exponential notation is adopted. 


p(x,y,t) = p (x,y) e 


io>t 


v-(x,t) = V (x) 

a a 


Dt 




where u is the frequency of oscillation and v the downwash on the airfoil surface. The 
boundary value problem is posed as follows. ^ 


V^p.o 

1 3e “’a , , 

- on the surface 

2 2 

p-»- 0 as (x + y )-►» 

p singular at the leading edge such that 
V (x,t) = v_(x,t) on the surface 

o 


(FI) 
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The pressure gradient condition ensures that the fiow has the correct curvature on the 
surface. It is then sufficient to satisfy the normal velocity condition at any one point, such 
as the leading edge. 

The complete solution can be separated into a singular part and a regular part 


p = p . + p 

^ *^sing *^reg 

The singular part, Pgjpg* which satisfies the condition " p |^ = 0 the surface, is singularj^ 

the leading edge. The regular part, P^eg’ continucy^ everywhere and satisfies “ ^ 

on the surface. For an airfoil at rest in steady flow, 0> snd the solution contains only a 

singular part. In terms of elliptic coordinates ( )» the singular solution can be shown to be 

(see ref. 17) 


° ^ ^h - fAos — 

where g is a constant to be determined. The elliptic coordinates are described by the 
relations 

X = b cosh n cos (j) 
y = b sinh r) sin (j) 

The regular solution can be obtained by separation of variables 

pU^ 

Substitution into the Laplace equation yields the ordinary differential equations 

H" - n^ H = 0 
2 

$" + n $ =0 
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for which the general solutions are 


H = Cje"’l + C2e""^ 


sin n 4> + cos n <j) 


For the solution to be periodic in <J> and antisymmetric about <p = 0, n the sine solutions 
above are admissible for (j) and n must be an integer. For the solution to die out at infinity, 
only the negative exponential solution can be used for H. The regular solution is then of the 
form 





a„ e”*' sin n 6 
n ^ 


(F3) 


Applying the pressure gradient condition, 


1 


9P. 


reg 


pU^ 3n 


n=l 





from which the coefficients are obtained as 



(F4) 
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where k = w b/U is the reduced frequency. 

It only remains to determine the constant g and this is done by applying the normal 
velocity condition 


1 3p I 

■pj SYolyo” 

-oo 


dto = Va(x,t) 


The integration is conducted along the linearized trajectory of the fluid particle, given by 


^o = ^ ^^o 


y 


o 


= 0 


Hence, 


(oto = wt+g(Xo-x) 


= p(x^,,yo>e 


io) t ik (x -x)/b 
e o 


Integrating up to the leading edge ^ = -1, 
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Substituting for the pressure, and letting 


f„(k) 5 ^ e-"" dn 

0 


the condition becomes 





The integral on the right is divergent and it is shown in Appendix H that 


Finite Part 



-ik cosh n 

€ 

(cosh n -1) 




which yields 



(F5) 


The solution for the pressure field is complete and associated quantities can be derived, e.g. 
Sectional lift: 


ii 

pU^b 


2ir 



(-sin (|) ) d (j> = -2 Tr(g + ^ a^) 


(F6) 
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The final expression for the surface pressure differential is 



^upper 


^ower 





n=l 



n=l 


a^ sin n <j) 


(F7) 


where x = cos <j). 

This formulation will now be applied to two specific cases of oscillatory downwash so 
that the results obtained may be compared with iavailabie results obtained using the velocity 
potential approach. 

(1) Harmonic pitching and heaving: 

If h and a denote the amplitudes bf the heaving and pitching motions and a the pitch axis 
iocation, 


Vx) 

U 



-Ik CR + (x-a) a] 


It can be shown that 


a^ = -2 ikci + km + a 5 ) 

_ k^g 
^2 - 4 


a^ = 0 (n>3) 


fl =-|Hj^2^k)-- 


-ik 

liT 
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-ik 



-f[«i 


( 2 ) 


(k) + 



+ 7T 



(k) 


y-W r- T 

— Q- = -a-ik Lh- (1 + a) aj 


It follows that 


g = C(k)[-o-ikR-ika( 5 -a)] +i^ 


H.^%) 

where C(k) = — — 


[2C(k){-5-il«*ik5(a.i)) 


- 3ik a 




(F8) 


This result can be shown to be identical with that calculated by the vortex approach (ref. 30, 
eq. 5-342). 
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(2) Stationary airfoil in a steady stream with a sinusoidal vertical gust: 
denotes the amplitude of the gust, then 


Va(x,t) 

U 


-^e-ikXeio)t 



0 


Hence the problem contains only a singular solution and it can be shown that 


- = 

^ ink 



pU^b “ ‘ ik[Hj^%) + iHQ^%)] 

The corresponding result from the velocity potential approach (ref. 30 eq. 5-376) is 

= 2 ir Vg [C(k) j 3Q(k) - i3^(k) j + i3j(k)] 

The two results appear to be different, but the second can be reduced as follows. 

C(k) [3Q<k) - i3j(k)] + i3^(k) 

= ~( 2 ) — ' — ( 7 ) — 


(F9) 


(FIO) 
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Using the recurrence relations for cross products of Bessel functions (ref. 31, 9.1.32 and 
9.1.34), the expression inside the square bracket can be shown to be equal to (2/ k). The 
result for the lift becomes 



which is identical with the result of the pressure method. 
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APPENDIX G 


HIGH ASPECT RATIO WING WITH OSCILLATORY DOWNWASH 


The coordinate system used for the straight, rectangular wing is shown in figure U. 
The linearized boundary value problem can be stated as 


V^P = 0 


2 2 2 

p-»- 0 as (x + y +2 )->oo 


* , Dv 

^ ~ on the wing surface 


(Gl) 


p->- 00 along the leading edge such that 


t 

-00 


on the wing 


A solution is required for wings of large aspect ratio, A. The aim of the analysis will be to 
obtain an approximate solution by neglecting all terms of order greater than 0 (A"^). 

The near field solution . - In the vicinity of the wing surface, excluding the tip regions, 
it is assumed that the characteristic length for spanwise variations is the span while the 
characteristic length for variations in the other two directions is the chord. The Laplace 
equation is then written as. 


2 2 2 
912- + 1 I 2 _.J_ rje. 
2 2 ~ 2 2 

<s) ^ 4f 


where, by the above assumption, the partial derivatives are of the same order. In the limit 
A -><» , the near field is a solution of the two-dimensional Laplace equation. For the large 
aspect ratio wing, the asymptotic expansion 
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P = Po ^ A Pi • • • 



shows that, neglecting terms of order greater than 0(A“b, the near field is still described by 
the two-dimensional Laplace equation, 


3 p 3 p 

*^near *^near 

2 + 2 
3 


with the conditions 


3P, 


near 


3y 



on the wing leading edge such that v(x,z t) = v (x,z,t). As shown in 

Appendix F, the^neral solution is 


1 

py2 ‘^'^ear 


= g(z,t) 


sin <}> 

coshTi + cos <J) 


(G2) 

+ ^ faj^(z,t) + b^(z,t) e'^’^ J sin n(f) 

n=l ^ 


the positive exponent being retained since there is no condition at infinity. The coefficients 
are given by 
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It can be shown that the far-field behavior of the near field solution is described by 



n=l 


where r and are cylindrical coordinates. 

The far field solution. - At large distances from the wing, the characteristic length can 
be taken to be the span for variations in all directions. The far field problem is p. singular 
along the line, x = y = 0 (-s s z < s), and antisymmetric relative to the plane y = 0. 


’ Piar = “ 


2 2 2 
P|^j. .►O as (x + y + z 


(G3) 


As shown in Appendix B, the far field solution to 0(A’b is that of a line of dipoles along the 
mid-chord of the wing. 


Pf 


ar 


pU^ 


1 

-1 0 


(G«f) 


where the starred variables are non-dimensional with respect to the semi-span, s. 


Lim 

* 

r o 


K,(qr*) = -V 

qr 
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The near field behavior of the far field solution is described by 


Urn 1 r f ( 5 *t)dc* r cos jq(c*-zM dq 


sin Y 
* 

2ttt 


f j(z ,t) 


The matching condition . - In order to construct a composite solution that tends 
correctly to the near and far field solutions in the respective regions, and varies smoothly 
from one to the other, it is required that 


V™ Pnear = V" '“‘ar 

I* 0 


By inspecting the limiting behavior of the solutions, it is seen that 


fj(z*,t) = + j aj(z*,t)] 


b (z ,t) = 0 for all n 
n 


The composite solution is written as 


P ~ Pnear Pfar ” Pcommon 
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where 


p = Lim p 

^common *^near 

X>*-00 


= Lim 
I.O 


Pfar 


The composite solution to 0(A”b is 


-2 = coshf4os » - E 


_2_ _ 
pu 


n=l 


1 00 

+ y I J* f aj( c *,t) J d C*J* qk^ (qr*) cos jq(c *-z*) j dq 


- j^g(z*,t) + 1 aj(z* t)] 


(G5) 


Since the problem pertains to oscillatory downwash, all quantities will be assumed to have 
complex exponential time dependence. The sectional lift is given by 


— ^ (z*) = -2 
pU^b 


1 

/ "^1 dx = -2Tr[i + i aj 

•{ pU I y = o 


which confirms the result of matching, viz that the dipole strength distribution is 
proportional to the sectional lift. The unknown function g (z*) is to be determined by 
applying the normal velocity boundary condition, 



y = o 
^o 


dto = v^(x,z,t) 
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where the integration is along the linearized trajectory 


yo = 0 

^o = ^ 

Integrating up to the midchord line, the condition becomes 



Substituting for the pressure gradient, letting 5 =--^ = Ar* and rearranging, the equation 
can be written as 


-2 

i-nk 


[ V (0,z ) « n - * 

n=l 


n=l tt72 


TT 


fi,(z*) 
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The left hand side, when multiplied by -2n, is seen to be the expr^gion for sectionai iift 
from two-dimensional theory (Appendix F),‘and will be denoted by-^^-y (z*). The first term 


on the right is the actual sectional lift, to the order of approximilior!^ considered here. 


VIZ 


0(A"b. The second term on the right can be put into a more convenient form, as shown in 
Appendix H. The resuit is the following integral equation for the sectional lift distribution. 


i^) £(z*) 

pu\ pU^b 



fe-iKSf A(.V) 

i e [<fFTP^ - 


( g -2 ) 

1 * *1 

|?-z |J 


d5 


X 


(G6) 


The integral term represents the effect of finite span to 0(A 


-1 


). In kernel form, 


fV) 

pU^b 


Kz*) _ 1 

pu\ ink TI^^^TTh^ 



— * (-^U[kA(cV)]d,* 

dC ^pUV 


where 


K(x) = k 



64 



Integrating by parts, 


K(t)=- 




- k T 


00 

X ItI 


- J 

00 

. f . 

d X 



dX 


= -ikN^(T) -kNg(x) 


(G7) 


The functions Ng and Np are the same as those discussed by Ashley, et al. in reference 32 
and can be expressed in terms of special functions. 


Nb(t> = i - i . i|[l,(|,|).Lj(|x|)]W 

-[■^K,(|t|) -i] 

N^(t) = Y + if + ln(2 |t|) 
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where I and K are the modified Bessel functions of the first and second kind and L the 
modifie9 Struve^Vunction. 

Reduction to the case k=0 . - In the limit k-»- 0, the following limiting values are valid. 


i iTk + 1Hq^%)] -► -2 

N3(t)- i-i 


N^,(t) ->y + ln(2 |t| ) + i| 


K(t)-> -f = - 


A( C -z ) 


The integral equation becomes 


iPiz) _ a(z*) _ _L f _d_ / I \ 
pU^b pU^b J d?*VU^b/ 



* * 


? -z 


which is the same as Prandtl's result for a high aspect ratio wing in steady flow. 

Comparison with Reissner*s result . - Reissner (ref. 19) represents the wing by a 
distribution of vorticity on the planform and in the wake, and then applies standard lifting- 
line approximations to the downwash integral. The basic result of this analysis is an integral 
equation for the spanwise variation of the reduced circulation. 


U 


-/ *v 

_ ^(z ) 
U 


3Q(k) -iJ.(k) r . p * * 


where ji=e*^ f and f is the bound circulation. Unlike steady flow, the lift distribution in 
unsteady flow is not completely determined by the bound circulation. However, the result of 
the pressure method can be compared numerically with Reissner's result. Numerical results 
for the latter, for the case of a uniformly pitching rectangular wing, are presented in 
reference 20, and the comparison is tabulated as follows. 
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Reduced 

frequency, 

k 


.333 
.333 
0.333 


0.667 

0.667 

0.667 


Spanwise 

Asymptotic lift 

coordinate. 

derivative. 

z* 

SjLsl 


pU^b 




Aspect ratio, A = 3 


Aspect ratio, A = 6 


Reissner's lift 
derivative, 

?■/ ri 

pU^b 


.0 

4.087 + 0.3551 

4.126 + 0.4961 

.4 

3.899 + 0.3621 

3.926 + 0.4541 

.8 

2.961 + 0.3241 

3.002 + 0.3141 

.0 

3.848 + 0.9931 

3.914 + 1.0721 

.4 

3.693 + 0.9861 

3.746 + 1.0821 

.8 

2.894 + 0.8351 

2.946 + 0.9481 

.0 

3.625 + 2.5371 

3.674 + 2.6201 

.4 

3.533 + 2.4901 

3.566 + 2.6021 

.8 

3.049 + 2.0911 

3.034 + 2.2581 


.0 

4.659 + 0.0831 

4.668 + 0.0861 

.4 

4.533 + 0.1291 

4.558 + 0.1561 

.8 

3.712 + 0.2441 

3.810 + 0.2781 

.0 

4.203 + 0.7921 

4.242 + 0.8021 

.4 

4.120 + 0.8331 

4.146 + 0.8601 

.8 

3.479 + 0.8721 

3.570 + 0.9361 

.0 

3.809 + 2.5081 

3.836 + 2.5161 

.4 

3.766 + 2.5171 

3.782 + 2.5441 

.8 

3.389 + 2.3531 

3.450 + 2.4521 


It is seen that the two sets of resuits are quite close, although it is difficult to associate the 
observed deviations between the two possible causes, viz the numerical calculation process 
and actual differences in the two equations. 

An unsteady lifting line, theory has also been derived by James (ref. 33), using the 
acceleration potential approach. This result should be directly comparable to the one derived 
here. However, it has been pointed out recently (ref. 34) that James' results are in error and 
hence this comparison was not carried out. 
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APPENDIX H 

miscellaneous integrals 


Evaluation of the integral (see Appendix F) . 


/ -ik cosh ri 
coshn-l 


As it stands, the integrai is divergent, but the divergent part can be isolated by writing 


f l-i'^cosh n^-lk) f ^-ik^ „ 

J cosh n “ 1 ^ y cosh n - 1 


= I, . I2 


where I 2 is divergent. The integrand in I. is indeterminate, of the 0/0 form at the lower 
limit, but application of L' Hospital's riile ^ows it to be finite. 


^ f (c-i)Vr-i 


Using 


dg — - 

( g-1) Vg -1 



-d Vx/|rf 
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00 

-ik J 
1 




Since 



1 


The first integral on the right can be expressed in terms of Hankel functions 

Ij = ^ [Hj^^^k) + iHQ^^^k) j + e"^*^ 


The finite part of can be extracted by recognizing it as the integral that would appear in 
the case of steady flow (ref. 17). This value of the finite part of I2 is -1, thus 


I = Ij + e“^*^(-l) 

= ^ [Hi^^\k) + iHQ^%)] 
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The integral with respect to u is to be interpreted as 


V 

f 


00 


o 




f (u) du = Urn / f (u) e’ ^ “ du 


<y. M I ^[cosff (c‘-z*)j dc 


( c *-z*) d [cos j Y ( j]j 


Making use of Z (+1) = 0, 


00 _ 00 
-ik? 


<Vi=-i jK(u)]cosjf(cV)jdu 

-1 ^ 0 ^ 0 


The inner integral can be evaluated as (ref. 31) 


f" K„(u) cos { !^ ( C V)j dd .= I [l * ( f ) " ( c *-zV ]■ 

“ [ I tH y i ftv| “ '*] 


1/2 


Ifs = 


.1 ^ 

J dc ‘ 

-1 



dc* 


rVinf ■ 

J 


A( g -Z ) (C -2 ) .j. 

C -z I ? -z 
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Proiate-sheroidal 
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= r cos X r _ s sinh ^ sin 0 

y = r sin X = s cosh cos e 
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z 
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Figure 3. - Cylindricai and Prolate-spheroidal coordinates. 
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Figure 5. - Normalized total blade lift versus azii 
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position for Case 1. 













































>- Measured 


® Asymptotic 



Spanwise location, r/R 

(k) Azimuth position, 300*^. (1 ) Azimuth position, 330°. 

Figure 7. - Concluded. 
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Sectional lift * span / thrust per blade 


Measured 


® Asymptotic 






span / 


Measured 


o Asymptotic 



Spanwise location, r/R 

(g) Azimuth position, 180®. (h) Azimuth position, 210®. 

Figure 8. - Continued. 
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Differential pressure * chord * span / thrust per blade 


Measured 


o Asymptotic 















Differential pressure * chord * span / thrust per blade 


Measured 


® Asymptotic 















Differential pressure * chord * span / thrust per blade 


Measured 




(b) Azimuth position, 45°. 



(c) Azimuth position, 90°. (d) Azimuth position, 135°. 

Figure 11. - Differential pressure versus chordwise location for CaseT 
with advance ratio of 0.29 at 95 percent span. 
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— Measured 
O Asymptotic 
A SVW 
^ Disk 


1 


i 


Azimuth angle, deg 
(a) Advance ratio, 0.06. 

Figure 12. - Normalized total blade lift versus azimuth position for Case 2 




Total blade lift / thrust per blade 



Azimuth angle, deg 
(b) Advance ratio, 0.13. 
Figure 12. - Continued. 







8 

Azimuth angle, deg 
(c) Advance ratio, 0.29 
Figure 12. - Concluded 




Sectional lift * span / thrust per blade 
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Sectional lift * span / thrust per blade 


Measured 


® Asymptotic 



Spanwise location, r/R 

(g) Azimuth position, 180°. (h) Azimuth position, 210°. 

Figure 13. - Continued. 
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Sectional lift * span / thrust per blade 


Measured 


o Asymptotic 





Sectional lift * span / thrust per blade 


Measured 




(g) Azimuth position, 180°. 


® Asymptotic 



WBSiK^ 


H^Hm 



0 .2 .4 .6 .8 1.0 

don, r/R 

(h) Azimuth position, 210°. 


Figure 14. - Continued. 
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Sectional lift * span / thrust per blade 


















Differential pressure * chord * span / thrust per blade 



(a) Azimuth position, 0°. 


o Asymptotic 




Figure 16. - Differential pressure versus chordwise location for Case 2 
with advance ratio of 0.29 at 75 percent span. 
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Differential pressure * chord * span / thrust per blade 


12 


Measured 


o Asymptotic 


9 [ 



-3 




(f) Azimuth position, 225°. 



0 .2 .4 .6 .8 1.0 

Chordwise location 

(g) Azimuth position, 270°. (h) Azimuth position, 315°. 

Figure 16. - Concluded. 
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Differential pressure * chord * span / thrust per blade 


Measured 


o 


Asymptotic 



Chordwise location 

(c) Azimuth position, 90®. (d) Azimuth position, 135®. 

Figure 17. -Differential pressure versus chordwise location for Case 2 with 
advance ratio of 0.29 at 85 percent span. 
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Differential pressure * chord * span / thrust per blade 


12 


Measured 


o Asymptotic 



Chordwise location 

(c) Azimuth position, 90°. (d) Azimuth position, 135°. 

Figure 18. - Differential pressure versus chordwise location for Case 2 
with advance ratio of 0.29 at 95 percent span. 
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Differential pressure * chord * span / thrust per blade 


Measured 


® Asymptotic 



Chordwise location 

(g) Azimuth position, 270°. (h) Azimuth position, 3i5°. 

Figure 18. - Concluded. 
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Total blade lift / thrust per blade 



Azimuth angle, deg 
(a) Advance ratio, 0.29. 

Figure 19. - Normalized total blade lift versus azimuth position for Case 3. 




Azimuth angle, deg 
(b) Advance ratio, 0.39. 
Figure 19. - Continued. 



Total blade lift / thrust per blade 



Azimuth angle, deg 
(c) Advance ratio, 0A5 
Figure 19. - Concluded. 





Total blade lift / thrust per blade 



Azimuth angle, deg 

Figure 21. - Computational effect on asymptotic total lift for Case 2 with advance ratio of 0.29. 


VO 


Total blade lift / thrust per blade 


N) 

O 


O 



Figure 22. - Computational effect on asymptotic total lift for Case 3 with advance ratio of 0.39. 
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